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, Abstract. In this paper we study the electromagnetic fields generated by a Killing 

,_H ' vector field in vacuum space-times (Papapetrou fields). The motivation of this work 

^ ^ ■ is to provide new tools for the resolution of Maxwell's equations as well as for the 

search, characterization, and study of exact solutions of Einstein's equations. The first 
part of this paper is devoted to an algebraic study in which we give an explicit and 

I . covariant procedure to construct the principal null directions of a Papapetrou field. In 

^ ' the second part, we focus on the main differential properties of the principal directions, 

■ studying when they are geodesic, and in that case we compute their associated optical 
, scalars. With this information we get the conditions that a principal direction of the 
' Papapetrou field must satisfy in order to be aligned with a multiple principal direction 
. of the Weyl tensor in the case of algebraically special vacuum space-times. Finally, we 
I illustrate this study using the Kerr, Kasner and pp waves space-times. 

"o ' 

^ ; PACS numbers: 04.20.-q, 04.40.Nr 

> ■ 

. 1. Introduction 

' In a relevant and historic paper, Papapetrou |jl|] pointed out that a Killing vector field 
(KVF), say can be considered as the vector potential generating an electromagnetic 
field with current = and satisfying a covariant version of the Lorentz gauge. 

As an immediate consequence, in vacuum space-times this electromagnetic field, which 
we will call here the Papapetrou field, satisfies the Maxwell equations in the absence of 
electromagnetic sources. 

The Papapetrou fields have been used in the search and study of exact solutions 
of Einstein's field equations (see e.g. 0, |^), and in the study of black holes under 
external electromagnetic fields 0. In particular, it turns out that the Kerr-Newman 
electromagnetic field is a Papapetrou field generated by the timelike KVF of the Kerr 
metric (for more information see references @, 

The applications and motivation of this paper concern with issues in which the 
Papapetrou field plays a central role. The first one is the resolution of Maxwell's 
equations in a curved space-time, which is a subject of special relevance for instance 
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in the study of electromagnetic perturbations of black holes (an exhaustive account is 
given in P). When write Maxwell equations in a null basis, using the Newman- Penrose 
formalism |^ (see also ||^, |^), we get a set of first-order partial differential equations 
for the complex components of the electromagnetic field, namely $a (^ = 0, 1, 2). The 
integrability conditions for the complex component $i (the conditional system) form in 
general a third-order differential system (with respect to the components $o aiid $2), 
however, as it was shown in if we choose the Newman-Penrose basis to be adapted to 
the principal null directions (or null eigendirections) of any particular regular solution 
of Maxwell's equations, the resulting conditional system is of second order. This is 
an important fact in order to integrate Maxwell's equations. For instance, as it was 
remarked in 0], this is just what happens with the Teukolsky-Press relations |1T0| , 
which were the starting point for the integration of Maxwell's equations in perturbed 
Kerr space-times (see and references therein). In [^, it is shown that the Teukolsky- 
Press relations, completed with an additional equation, constitute a conditional system 
of second order for Maxwell's equations (in the variables $0 and $2)- This is due to 
the fact that they were written using a Newman-Penrose basis adapted to the principal 
directions of the Kerr space-time, which are also the principal directions of a regular 
electromagnetic field (see also P]). These facts emphasize the importance of knowing a 
particular regular solution of Maxwell's equations. In general space-times we do not have 
a method to construct such a solution, but in the case of vacuum space-times possessing 
a KVF, a particular solution is given by the Papapetrou field. Hence, part of this paper 
is devoted to make an exhaustive algebraic study of these electromagnetic fields, giving 
the procedure to construct a Newman-Penrose basis in which the integrability conditions 
of the Maxwell equations are directly second-order partial differential equations. For 
the sake of completeness we will also treat the case of singular Papapetrou fields, which 
appear in relevant space-times like the well-known pp waves. 

On the other hand, Papapetrou fields provide us a link between the Killing 
symmetries and the algebraic structure of the space-time, which is a subject treated 
only occasionally in the literature (see Chapter 33). This relationship can be found 
by studying the alignment of the principal directions of the Papapetrou field with those 
of the gravitational field, i.e., of the Riemann tensor (which in vacuum reduces to the 



Weyl tensor). In this sense, a remarkable example is the Kerr metric [jT2|. In [T^, ^ 



it is shown that the principal directions of the Papapetrou field associated with the 
timelike KVF coincide with the two repeated principal directions of the space-time (it 
is Petrov type D). This link between symmetries and algebraic structure provides a 
powerful tool for the study and search of vacuum gravitational fields. As an example, a 
new characterization of the Kerr metric can be found. In a recent work by Mars |]14[, it 
is shown that the Kerr metric is the only stationary asymptotically-fiat vacuum space- 
time in which the aforementioned alignment is given. This characterization of the Kerr 
metric is obtained everywhere, not only where the Killing is timelike. Here, we will see 
an alternative and simple way of characterizing the Kerr metric by using our formalism 
and the connection with the eigenray formalism of Perjes ||T5 |. 
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These facts motivate the study of the differential properties of the principal 
directions of the Papapetrou field. In particular, the study of under which conditions 
they are geodesic and shear-free will tell us, via the well-known Goldberg-Sachs 
theorem |]T6[, when a principal direction of the Papapetrou field is aligned with a 



multiple principal direction of an algebraically special vacuum gravitational field. As we 
will see in the development of the paper, these properties, and other properties related 
with the optical scalars of the principal directions, can be determined only in terms of 
the principal direction itself and the Ernst potential associated with the KVF under 
consideration. This result is relevant in the sense that it provides a way of making 
interesting Ansatze for the search of new solutions of Einstein's equations which, in the 
light of the example of the Kerr metric and other examples that will be given in this 
paper, can be physically well motivated. 

The plan of this paper is the following: in Section H, we give some general properties 
of antisymmetric tensors (2-forms), which will be needed for the development of the next 
sections. In Section ||, we introduce the Papapetrou field associated with a KVF. Then, 
we make an exhaustive study of its algebraic structure, for which we shall distinguish 
between regular and singular fields. In the regular case, we will introduce a complex 
1-form, proportional to the differential of the Ernst potential, which will play a central 
role in our study. Then, we will determine explicitly (in terms of quantities defined only 
from the KVF and the metric) the eigenvalues, the principal null directions, and the 
orthogonal 2-planes to these directions (the Maxwellian structure). In the singular case 
we determine the only principal null direction. In Section |^, we use this algebraic study 
to analyze the main differential properties of the principal direction(s). We will give the 
conditions for a principal direction to be geodesic, and for that case we will give the 
expression for the optical scalars, which can be written in terms of the complex 1-form 
mentioned above and the principal direction. An illustration of how this study works 
is given in Section |^, where some examples are discussed, in particular the case of the 
Kerr metric [|l^ is treated with great detail. In Section ^ we discuss the main results 



and consequences of this work, as well as possible extensions. Finally, in Appendix A we 



summarize the known results about the case of a null Killing vector, and in [Appendix B 



we give the connection with the eigenray formalism of Perjes [0. Through this paper 



we will follow the notation and conventions of unless otherwise stated. 
2. Algebraic structure of a 2-form 

In this section we review briefly the algebraic structure of a 2-form F (an antisymmetric 
tensor), describing an explicit and covariant procedure to construct the null principal 
direction(s) (eigendirections) . First, we deal with the regular (non-singular) case, that 
is, when at least one of the two invariants is not zero: 

F'^'Fab ^ , Fab = Fab + i* Fab , *Fab = ^VabcdF''' , 
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where * and ~ denote the dual and self-dual operations respectively, and rjabcd are the 
components of the volume 4-form. At the end of this section we will extend the procedure 
to the singular case. 

A regular 2-form F can be completely characterized in terms of its eigenvalues 
and null principal directions {k,£). In fact, it can be written as follows 

F = aG- (3 *G, G = kA£, fcX = -1 • (1) 

Note that G is a 2-form of rank 2 [det(G) = 0]. Moreover, the corresponding dual and 
self- dual 2-forms *F and F are given by 

*F = (3G + a *G, F = {a + i(3)G, (2) 

from where we find the following relation between the invariants of F and its eigenvalues 

F'^'Fab = -4(a + 1(3)'' . (3) 

Given a 2-form F, this relation allows us to find the eigenvalues (q;,/3), and from (|^) we 
can find the singular 2-form G. On the other hand, for each 2-form F we can construct 
the following symmetric tensor 

Tafe = ^(F/Ffe, + *F/*F,,) . (4) 

When F represents an electromagnetic field, this tensor is the energy-momentum tensor. 
In the case of a regular 2-form, using (|l],0) we find an alternative expression for Tab 

Tab = (a^ + P') k{ah) + rn^arrib) , (5) 

where m is a complex vector orthogonal to the 2-planes generated by the principal 
directions k and ^, and such that {fc,£,m,m} is a Newman-Penrose basis. Therefore, 
this basis is adapted to the Maxwellian structure of the 2-form F. We can write the dual 
of G in terms of m and its complex conjugate m as *G = imAffi. Note that expression 
(^) shows the 2+2 decomposition of the energy-momentum tensor (equivalently of the 
Ricci tensor) of a regular 2-form, or in other words, that the Segre type is [(1 1) (1 , 1)] 
(see for instance 

The problem of determining the principal directions of a 2-form has been treated 
by Coll and Ferrando fl^ • To that end they introduced the so-called concomitant of a 
2-form 

J^ab = aFab + P* Fab-Tab + ^ + P^)9ab ■ (6) 

Using this tensor, the principal null directions are given by 

Ka = J'abU' , La = J'baU' , (7) 

where U is an arbitrary timelike vector field. We can understand this property by 
rewriting the tensor in the following form 

= («' + /9') {Gab + Ga'G,,) = -2{a' + P')iakb , (8) 



where we have used (^^- As we can see from this expression, K and L are parallel to 
the null vector fields k and i respectively [see expressions (|l|)]. It is important to point 
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out that this method provides a covariant way of finding the principal directions of a 
2-form. 

On the other hand, in order to determine the 2-planes orthogonal to the principal 
directions, we can introduce another tensor 

= i{(3Fab - a * Fab) +Tab + l (a' + P^)gab . (9) 

From this complex tensor we can find a complex vector m describing the orthogonal 
2-planes. To that end, we have to contract (P) with any spacelike vector field linearly 
independent from k and £, and normalize the result {m°'fha = 1). We can show this fact 
by rewriting JF^ in terms of the Newman- Penrose basis {k,£, m, ffi} 

^ab = -(«' + * Gab + *Ga' * G,b) = -2(a' + (3^)mam^ , 
Note that m is fixed up to a factor e**", where C is an arbitrary real scalar. 

To sum up, we can construct the Maxwellian structure associated with a regular 
2-form from tensors (|^) and (^. 

In the singular case, the self-dual 2-form F can always be written as follows (see P]) 

Fab = 4:^k[amb] , k^rua = , 

where m is some normalized null complex vector field {m"'ma = and m'^fha = 1) and 
$ is a complex scalar field. As is well-known. A; is a multiple eigenvector with zero 
eigenvalue. In addition, the energy-momentum tensor @ now takes the following form 

Tab = 2<l>Makb , 

that is, the Segre type is [(1 1 , 2)]. We can use the previous procedure to find explicitly 
the principal direction only by considering the fact that the eigenvalues are identically 
zero {a = (3 = 0). Therefore, the tensor J^ab is simply —Tab, showing that there is only 
one principal direction [K and L coincide), which is determined by 

Ka = TabU\ 

where U"" is an arbitrary timelike vector field. In that case, the tensor JF^ does not 
provide any information since there is not a 2-1-2 structure, which is clear from the 
Segre type of the energy-momentum tensor. 



3. The Papapetrou field 



Let {V4,g) be an arbitrary vacuum (Rab = 0) space-time (see [|T8l for details) endowed 
with a one- dimensional isometry group generated by a non-null KVF ^ 

^a;b + 6;a = . (10) 

That is, we will assume that $ has a fixed character (timelike or spacelike) in an open 
domain of the space-time[[]]. In this domain, the electromagnetic field, Papapetrou field, 
generated by $, is given by 

Fab = {d$,)ab = Cb;a — Ca;f> = 2^fe;a ) (H) 

II The case of a null KVF is of little interest since all the vacuum space-times admitting a null KVF 
are known, they can be found in js) (section 21.4; see also Ref. jl^ and Appendix A ). 
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where we have used the Kilhng equations (|TOD and d denotes exterior differentiation. 
Note that F is defined up to a multiphcative constant. Moreover, taking into account 
the following relation for the second covariant derivatives of a KVF 

^a;bc = Rabcdi'^ , (12) 

which comes from the Ricci identities, we can show that this 2-form Fab satisfies the 
Maxwell equations in the absence of electromagnetic charge and current distributions 

F[ab;c]=0, F'^% = 0. (13) 

In addition, from the definition of the Papapetrou field (|lTl), we can see that the KVF 
$, plays the role of the electromagnetic vector potential, and from the Killing equations 
(p!0|), that it satisfies the covariant Lorentz condition ^"^.^ = 0. 

Some important quantities associated with a KVF are: the norm 

which we have assumed to be different from zero, the orthogonal projector 

hab = Qab — J^iaib , hab^^ = , 

the gradient of the norm 
and the twist 

The twist has three independent components which are also contained in the following 
2-form 

Wab = ha'hb%c;d] = K'hb%d , Wabe = , (14) 

usually called the rotation of The rotation is related to the twist through the 
expressions 

OOa = Vabcd^W'''' , Wab = ^Vabcd^Oo" • 

In terms of these quantities, the Papapetrou field can be written in the following 

way 

2 2 

Fab = ^ {C[aA] - NWab) , *Fab = — {^[a^^b] + NMab) , (15) 

where the definition of Mab is 

Mab = ^VabcdCi^"^ ■ 

On the other hand, it is well-known that in vacuum space-times the twist 1-form 
satisfies du; = (see for instance), and therefore there is locally a function Q so that 
= dfl. Then, we can introduce the Ernst potential (see |l20| , §) associated with ^ 

£ = -N + in. 
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In terms of this potential, the self-dual Papapetrou field can be cast into the following 
form 

F = ^ A e)~ e = d£ = + iu: {de = 0) , (16) 

where we have defined e as the exterior derivative of the Ernst potential S. In vacuum, 
this 1-form satisfies the following equation 

e% + iV-ie% = 0, (17) 

which comes from Einstein's equations. Moreover, from (|T6|) we realize that in the case 
of a timelike KVF, tj) and u} are respectively proportional to the electric and magnetic 
fields measured by observers with 4-velocity u = (—N)^^^"^^. 

We begin the process of construction of the principal directions of the Papapetrou 
field by the calculation of the eigenvalues a and (3. Using ([T^) or (|T^) we get 

And comparing this expression with (H) we find 



s / rr. „ s 



a = + -y/x^ + 2/2, p = ^ — ^jx+ + 

where x and y 

— N — ' y^^^ (19) 

are the two invariants of the Papapetrou field F, and s is an arbitrary sign (s^ = 1). 
This sign reflects the invariance of F under the change {k,i, a, (3) (£, fc, —a, —(3). 

From (|TH]), we find some interesting relationships between the eigenvalues and the 
invariants 

- (3"^ = -x^ 2a(3 = -y , + (3^ = .Jx^ + y^ . (20) 

These relations together with (|19|) imply 

e"ea = N{x + iy) = -N{a + i^f . (21) 

On the other hand, from (^ we have the following expression from the energy- 
momentum tensor of the Papapetrou field (it is important to remark that this 
electromagnetic field is a test electromagnetic field) 

lN{ri^, + u''uc)gab] , (22) 



2 



where 



2Wabi^' = -2M,b^' = -Lrjabcaeruj" 



In this situation, we are ready to construct the tensors J-'ab and J-"^ (|^J^)) and 
through (1^ we can find explicit expressions for null vector fields tangent to the principal 
directions of the Papapetrou field. However, in order to carry out these calculations. 
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it is necessary to consider separate cases. First of all, we make a division into two 
groups depending on whether the Papapetrou field is regular {a + i(3 ^ 0) or singular 
{a + if] = 0). In the first case, we will determine completely the characteristic 2+2 
structure of a regular 2-form, the Maxwellian structure. In the second case we will give 
the only one principal null direction. 



3.1. Regular Papapetrou fields 

We begin the study of regular Papapetrou fields by introducing a complex 1-form which 
will be of crucial importance for the purposes of this work. The definition of this 1-form 
is 

E = — ^ = — ^(il) + iuj), EAdE = 0. (23) 

As we can see, it is proportional to the exterior derivative of the Ernst potential and 
hence, it is integrable. Moreover, from ([T6|) we find the following useful relationship 
between E and the self-dual 2-form G [see Eq. (^] 

G = ^{iKEr ^ Ea = -Gate- 

On the other hand, E satisfies the following equation 

a + tp 

which follows from equation (jl^). On the other hand, when it is an exact 1-form, 
dE = 0, either a + ip is constant or it is a function only of the Ernst potential. 

The importance of this 1-form is twofold, first it will allow us to classify the different 
subcases that need a separate treatment, and second, it will serve to write our results 
in a compact form. Now, let us split E into its real and imaginary parts 

E = En + iEi . (24) 

From the definition (^) we can find the relations between {Eji,Ei) and 

««=^("* + ^' = ^(-W + -). (25) 

and the converse relations 

if) = aEu - [3Ei , uj = PEr + aEi . 



From (|23|j2^ , |2l|) we find the following important relation 

E^E, = E^^E^, - E^E,, + 2zE^^E,, = -N , (26) 
from where we deduce that Epi and Ei are orthogonal 

E^REia = , (27) 

and also, a relation between their norms and the Killing norm (note that E has a real 
norm) 

- E'}E,, = -N . (28) 
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From these properties we can derive some consequences concerning the character 
of the vector fields Eji, Ej, (p}. First, when ^ is timehke (A^ < 0), is obvious that 
Er and Ej are spacehke, and taking into account that 

r<Pa = ^ [irc^af - iri^a)ico'co,)] = -j^{a' + P'f{E'i,E^ME:,) , (29) 

it is clear that (p is also spacelike. On the contrary, when ^ is spacelike (A^ > 0), from 
(p8|) we have EjEj^ = N + E'j^E^ and then, Ej cannot be timelike, otherwise Eji would 
be also timelike, in contradiction with (pTl). Moreover, Ej cannot be a null vector either, 
otherwise taking into account (0), which implies that E'^E^ = —N < 0, and (]27|), Eji 



would be also a null vector and using again (|28D, ^ should be also a null vector, which 
contradicts our initial assumptions. Therefore, we have shown that Ei must be always 
spacelike (unless it vanishes). In contrast, Eji and can be timelike, or spacelike, or 
null. Furthermore, when En is a null vector, ^ is spacelike and 4> is also a null vector 
parallel to Eji. 

Now, let us introduce some vectors which are adapted to 2+2 decomposition of the 
energy-momentum tensor of an electromagnetic field. These vectors are En, Ej and 

P ^ 1 [{E'i,E^,)i + r>] , Q = 1 mE,,)i + D] , (30) 

where is a 1-form defined by = *{$,AEjiAEi). In local coordinates the components 
of D are 

The scalar products between these 1-forms are given by 

P'^Pa = -E^nE^ , Q^Qa = E^Ej, , (31) 

P-'Ena = P^Eja = P^Qa = Q'^ERa = Q^^Ej^ = , 

and hence, they form in general an orthogonal basis (unless E'^E^^ = or Ei = 0). 
The vectors P and En lie in the 2-planes generated by the principal directions. This 
fact can be seen from the following properties of JF^ 

T^,p' = :FtA = , :f^,p'' = ^f^.e'^r = o • (32) 

On the other hand, Q and Ej determine in general (when Ej ^ 0) the orthogonal 
2-planes to the principal directions. This fact follows from the following relations 

^ abQ^ = ^ ibE\ = , JF ^^Q" = JF a^i?" = . 

Moreover, since Ej is always spacelike and taking into account (0), we deduce that 
both Ei and Q are spacelike. 

In what follows we find explicit expressions for the principal directions of a regular 
Papapetrou field and determine the 2-planes orthogonal to these directions. Following 
the previous discussion we will distinguish between two subcases, depending on whether 
or not the norm of En vanishes. 
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3.1.1. Generic case E^j^E^^ ^ 0; This case represents the generic situation. From 
we deduce that either P or Eji must be timehke and hence, in order to obtain the 
principal directions we can use these vectors. However, in this case we can write in 
terms of the basis {P, En, Q, Ej} 

Tab = rnr 7-1 {Epui + Pa) {E Jib — Pb) , 

from where it follows that the principal directions of the Papapetrou field are determined 
by the following null vector fields 

K = P-Eu, L = P + Er. (33) 

These expressions also show the fact that in general P and En determine the 2-planes 
generated by the principal directions, as it has been pointed out before [see expressions 
(^)]. In order to determine the orthogonal 2-planes we have to use the tensor (j^), which 
in this generic case can be cast into the following form 

= ^^^{Eja - lQa){Eib + iQb) ■ 

EjEj^ 

Then, the orthogonal 2-planes are determined by Q and Ej. Moreover, as is clear, this 
expression requires EfEj^ ^ 0, equivalent to Ej 7^ 0. In fact, the case Ej = needs a 
separate study, which is given below. 

Finally, we want to stress the fact that the 2+2 structure is also manifested in the 
form of the energy-momentum tensor. In our case, we can see this feature by writing 
Tab using the basis {P, Eji, Q, Ej} 

Tab = ^(«' + P^) (P-Pb - ERaEnb) + (QaQb + EjaElb)] • 

• Particular case Ej = 0: Now we have + aoj = and then if) and oj are linearly 
dependent, hence D = 0. Moreover, this implies P = and E = En, therefore 
we can use the previous expressions (|33D . Then, the principal null directions of the 
Papapetrou field are given by 

K = ^ + Er, L = i-En. (34) 

On the other hand, the relations EjEj^ = and E'^Ej^^ = —N [the last one comes from 
(H)] imply 

^|J-^l;^ + uJ-uJa = -N{a^ + P^), 

where we have used (|2B|). Combining the last equation with expressions ( pi?| , pUD we 
arrive at the following results 

r^a = -No" , u'^UJa = -NP^ . 

The peculiarity of this subcase is related with the structure of the orthogonal 2-planes. 
Now, although we can determine the 2-planes spanned by the principal null directions, 
we cannot determine explicitly the orthogonal 2-planes, since in this case we have 
Ej = Q = [see equations (|30|) ]. However, using the results of section 2 we can 
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construct the orthogonal 2-planes. We only need to contract JF^ (^), which is now given 
by 

^ab = {iriabcdP^Ej, + PaPb - EnaEm - Ng^b} , (35) 

with any spacelike vector field independent of ^ and Er. Then, after normalization 
we find m, which is defined up to a factor e**^, where C is an arbitrary real scalar. 
Alternatively, given two vector fields, linearly independent and orthogonal to ^ and Er, 
we can construct from them the orthogonal 2-planes to the principal directions, and 
therefore we can find a vector m. 



3.1.2. Particular case E'^Ej^ = 0: Now, let us see what happens when Er is a null 
vector 

E%Ej,, = =^ E^Ej, = N. (36) 

Following the discussion given above on the norm of Er and Ej, it turns out that this 
case can only take place when ^ is spacelike (A^ > 0). 

On the other hand, using (|25|) , equations (^) and (|36|) imply 

Combining this result with (p!9| , |20|) we get expressions for the norms of i/^ and a; 



then, both are spacelike. Moreover, from ([T9|j20|) we can also deduce the following result 

This is an important relation because it means that either i/j and lc are parallel or one 
of them is a null vector orthogonal to the other one. The first possibility implies that 
Er = 0, which will be treated later as an special subcase. The second possibility leads 
to one of the following two excluding situations: the first one is given by 

r^, = ij''ua = P = 0, Er=-iP, Ej = -LJ, (37) 

a a 

and the second one by 

oo-u, = r^J, = a = , ^R=^'^^ Ei = ~if). (38) 



In both situations ( |37| , |38D , equation (^) implies that is a null vector field parallel to 
Er, and obviously, the same happens with D. 

With regard to the principal directions, in this case we cannot determine completely 
the 2-planes spanned by the principal null directions, since P and Er are not linearly 
independent {P = N^^D is also a null vector parallel to Er). From (|33D we can see 
that Er is one of the principal directions, but we cannot determine the other one from 
these expressions. This is due to the fact that we have not a timelike vector available. 
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Therefore, it is convenient to follow the procedure explained in Sec. 0. First, we 
have to construct the tensor J^ab (H)- In this case, we can find J-'ab through the expression 
(§) and taking into account that now we can write 

Gab=^VabcdE'iQ''. (39) 

Finally, using an arbitrary timelike vector and expressions (0) we will find the remaining 
principal direction of the Papapetrou field. 

We finish this case with the study of the following remaining subcase: 

• Particular subcase En = 0: The vanishing of En = is equivalent to the relation 
at/) + jSu: = 0, which tells us that t/> and a; are linearly dependent. The difference that 
appears between this subcase and the general case E'^Ej^ = is that now we cannot 
determine explicitly any of the two principal directions, since now we have En = P = 
[see (|33|)]. Therefore, like in the previous case, we have to find a timelike vector field 
and use the tensor J-'ab (^, which in this subcase has the following form 

^ab = [^labcdQ^Ej + QaQb + EjaEif, — A^fi'afc} ■ 

Finally, the principal null directions are given by expressions (|^). However, if we have 
a timelike vector orthogonal to Ej and Q, say U {Lf^Eja = U^-Qa = 0), the principal 
directions can be found through the following expressions 

Ka = VabcdU'Q'^Ef - NUa , K = VabcdU'Q'Ef + NU^ . 



3.2. Singular Papapetrou fields 

As we have said before, this case is characterized by the vanishing of the complex 
function a + ijS, which is equivalent to the vanishing of the invariants of the Papapetrou 
field X and y (p!9|): 

t/-^, = u'^Ua , = . (40) 



Moreover, the exterior derivative of the Ernst potential e (16) is a null 1-form, e^e^ = 0, 
and the tensor J-'ab becomes —Tab, that is, it is symmetric, in agreement with the fact 
that for singular 2-forms there is only one principal direction. It is obvious that in this 
case we cannot introduce the 1-form E, but by virtue of (|40|), the exterior derivative of 
the Ernst potential e is a complex 1-form satisfying similar properties that E, in the 
sense that the real and imaginary parts of e are also orthogonal. 

In order to determine the principal direction of a singular Papapetrou field two 
subcases must be distinguished: (i) When t/? and u: are two spacelike 1-forms which are 
mutually orthogonal, (ii) The rest of possibilities, which are characterized by the fact 
that the 1-form e, which is a null complex 1-form, is proportional to a real null 1-form 
k. 
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3.2.1. Subcase (i): As we have said before, we consider t/? and uj to be spacelike and 
orthogonal. They cannot be timelike, since (^OD would imply that they should be null 
1-forms, which corresponds to the subcase (ii). On the other hand, either ^ or the 1-form 
(f) must be timelike. Then, we can determine the principal null direction by contracting 
the energy-momentum tensor with one of them. We find that this principal direction is 
given by 

k = r^,i-N4>. (41) 

3.2.2. Subcase (ii): This case, characterized by 

e = Afc , (42) 

where A is a complex function, includes some particular cases but in all of them if) and 
(jj are null or vanishing 1-forms. As we can see from expression (plSf ), the null vector 
k in (|4^) determines the principal null direction of the Papapetrou field. Finally, it is 
important to point out that in this case the principal null direction is orthogonal to the 
KVF ^, k'^ia = 0. 



4. Some properties of the principal directions of the Papapetrou fields 

We devote this section to study some properties and features of the principal null 
directions of the Papapetrou field in vacuum space-times. More specifically, we study 
under which conditions these null principal directions are geodesic, and in that case we 
compute the optical scalars associated with them. Here, we will also distinguish between 
the regular and singular cases. 



4.1. Regular Papapetrou fields 

First of all, we are going to study when a principal direction of the Papapetrou field is 
geodesic. We can do that by using the explicit expressions for the principal directions 
given in the previous section, or their properties in the case En = 0, since in this case 
we do not know them explicitly. Then, the condition of geodesicity for a principal null 
direction is equivalent to say that the components of the vector k^k°'-b on the orthogonal 
2-planes must vanish. Then, contracting this vector with the explicit expressions for 
the generators of the orthogonal 2-planes (excepting in the case = 0, since in this 
case we do not have explicit expressions for them) and making some straightforward 
calculations, we can establish the conditions that a principal direction k must satisfy 
in order to be geodesic. The main point in these calculations is the use of the following 
relation between P, Q and ^ 

Q-P = ^, (43) 

which follows from (^) . This expression allows us to transform terms with Q into terms 
with P [and hence, into terms with the principal directions through (^)] and In this 
way, we can arrive at the following theorem: 
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Theorem 1: Let (V4, g) be a vacuum space-time and let ^ be a KVF. Then, a principal 
direction k of the Papapetrou field associated with ^ is geodesic if and only if the 
following condition holds 

k'^k^Ea-b = when Ej , 
m°-{a + i(3)^a = when = . 

Here, some remarks are in order. First, m is an arbitrary complex 1-form lying 
on the orthogonal 2-planes and such that m A m 7^ 0. We can find it by means of 
the tensor JF^ (|9|j35|). Second, when En is a null 1-form we have seen that it is also 
one of the principal directions, and from expressions (|37|j38| ) we deduce that Eji is in 
addition integrable, and therefore, it is geodesic. Finally, it is important to note that in 
the case Ej 7^ 0, the information on the geodesicity of a principal direction is encoded 
in a complex scalar depending only on the principal direction itself and the 1-form E. 
In the case Ej = 0, it depends on the gradient of a + i(3 on the orthogonal 2-planes. 

When one of the principal directions is geodesic, we can study its optical scalars: 
expansion, shear and rotation (they are only well defined in the case of geodesic null 
vector fields). We begin with the shear of the geodesies tangent to a given principal 
direction. As is well known, if we choose the tangent vector field k to be affinely 
parametrized, then the shear scalar coincides with the quantity aa, where a is one 
of the Newman- Penrose spin coefficients (see [0, §])• But in general, we can take a 
Newman-Penrose basis associated with k, namely {k,£,m,ffi} , and then to compute 
the shear tensor Sat 

5. . HM.> - Ih.H-K, . ff. = 2™,„*, ,H'^ = 2, (44) 

where Hab is the orthogonal projector to the 2-planes spanned by {fc,^}. Then, taking 
into account that the shear tensor Sab is spacelike, the shear scalar S"^ = \S'^^Sab 
vanishes if and only if Sab does, and it coincides with aa when k is chosen to be affinely 
parametrized. 



For our purposes it is better to compute the shear tensor (0). In the case Ej 7^ we 

can take m to be parallel to Ej + iQ (with a proportional factor such that rrf'Tha = 1), 

whereas in the second case it can be found through the tensor (p|,|35D. Then, after 
some long but straightforward calculations and using here also the relation (|43|), we 
arrive to the following result 

r ^^k-E'E[a;b] iiEj^O, 
S^ = ZZ , where Z = I 

\ !ffi.ni''rn'ERa:b = if Ei = . 

Then, using this result we can enounce the next theorem: 

Theorem 2: Let (V4, (7) be a vacuum space-time and let ^ be a KVF. Then, a geodesic 
principal direction k of the Papapetrou field associated with ^ is shear-free if and only 
if the following condition holds 

k''E''E[a;b] = when Ei ^ , 
rrf-m^ERa-b = when Ej = . 



On the Papapetrou field in vacuum 



15 



Here it is important to point out that the condition for the case Ej ^ involves 
only antisymmetric covariant derivatives and therefore, this calculation does not need 
the use of connection, we can make it by using only partial differentiation. In fact, it is 
possible to show that this condition can be written as follows 

\k^Eh)E^ + A^r] {a + i(3)^a = 

that is, in terms of derivatives of the eigenvalues of the Papapetrou field. Finally, taking 
into account the Goldberg-Sachs theorem |]16|, which states that a vacuum space-time is 



algebraically special if and only if it contains a shearfree geodesic congruence, the next 
corollary follows immediately 

Corollary: Let (V4,(?) be a vacuum space-time and let ^ be a KVF. If one of the 
principal directions of the Papapetrou field associated with ^, namely k, satisfies the 
following conditions 

k''k''Ea-b = and k'^E^E^a-fi] = when Ej^O, 

'm°'{a + i(3)^a = and m'^m^Ejia;b = when Ej = , 
the space-time is algebraically special and fc is a principal direction of the Weyl tensor. 

Now, let us study the expansion and rotation in the case of geodesic principal 
directions of the Papapetrou field. These quantities are given by 

^ = ^H^'^ka-b , Ha''Hb'^k[c-d] = -2iwm[arhb] , 

where ^? and zu are the expansion and rotation scalars respectively. In order to compute 
them it is more appropriate to consider the complex scalar 

p= -{{} + iw) = -rrfm^ka-b, (45) 

usually called the complex divergence scalar, which is another spin coefficient in the 
Newman-Penrose formalism. The result of this calculation is 

(k-E'^E^a-M + k-ER^E,,b) i^Ei^Q, 

[{k^ER,)m''m^ER(^a;b) + {k'QWabm'^m^] if = . 

In the case Ej = we can identify immediately the expansion and shear scalars, they 
are 

^ = ^^^rn'^rn'ERi^a.b) , vj = -z^Wabm'^m' . (46) 

Here, it is interesting to note the direct relationship between the rotation Wab of $ 
and the rotation of the principal direction. With regard to the case Ei ^ 0, it is 
remarkable that when the geodesic principal direction k is also shear-free, the complex 
divergence is simply 

P = ^k^En'E^.b , 
and hence, the expansion and rotation scalars are given by 

^ = -;^k^En'Eju;b , w = -^k'^En'Eja.b ■ 



P 
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Finally, it is important to remark the role that the 1-form E plays also in the differential 
properties of the principal null directions of a regular Papapetrou field, since as we can 
see in all the expressions of this section, the optical scalars as well as the geodesicity 
condition are given only in terms of scalars formed from the principal direction k and 
the covariant derivative of E. 



4-2. Singular Papapetrou fields 

The general situation in the singular case can be summarized as follows: we have a 
vacuum space-time containing a singular 2-form (the Papapetrou field), which is a 
solution of the Maxwell equations (|13D- Therefore, the Mariot-Robinson theorem [^1], ^ 



(see also tells us that this is equivalent to say that there is a geodesic and shear-free 
congruence of null curves. In addition, the principal null direction of the Papapetrou 
field will be tangent to these curves. Finally, the Mariot-Robinson theorem tell us that 
the space-time is algebraically special and then, the following theorem follows 

Theorem 3: Let (V4, g) be a vacuum space-time and let ^ be a KVF. If the gradient of 
the Ernst potential e is a complex null 1-form, then the space-time will be algebraically 
special and the principal null direction of the 

singular Papapetrou field associated with ^ will be geodesic and shear-free. 

Taking into account that we are dealing with vacuum space-times, this theorem 
tells us that if the gradient of any KVF in a vacuum space-time is a null 1-form, then 
the space-time must be algebraically special and a multiple direction of the Weyl tensor 
coincides with the null principal direction of the corresponding Papapetrou field. 

On the other hand, for the subcase (i), the complex divergence p ( ^5| ) is given by 
the following expression 

In the subcase (ii), using the equation ([17|) , we can see that p = 0. 



5. Some examples 



In this section we apply our previous study to some examples. In particular, we study 
the Papapetrou field associated with KVFs in the following vacuum space-times: (i) The 
Kerr metric, (ii) The Kasner metrics, (iii) A subclass of the plane-fronted gravitational 
waves. In these examples we compute the principal direction(s) of the Papapetrou field 
and also discuss some of their algebraic and differential properties. 



In the first example we consider the well-known Kerr metrics [|T^. They constitute 
a family of stationary axisymmetric vacuum solutions of Einstein's field equations. In 
Boyer-Lindquist (BL) coordinates {t,r,9,Lp} (see |2^, |]), the line-element takes the 
form 



A 



{dt — asin^ 9dipY + 



[r'^ + o?)dLp — adt 



+ ^dr^ + p'dO^ 



(47) 
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where A and p are given by A = r^ — 2Mr + a^ and = r^ + a^ cos^ 6 respectively. Here, 
we must say that in the case < M^, the BL coordinates are vahd in the asymptotically 
flat regions (r+ < r < oo) or type I regions, in the type II regions (r_ < r < r^), which 
contain closed trapped surfaces, and in the asymptotically flat regions containing the 
ring singularity = {—oo < r < r_) or type III regions {r± = M ± (M^ — a^)^/^; 
see for more details). The only singularity of this metric is the ring singularity 
p^ = 0. The event horizons (A = <S=^ r = r±) and the axis of symmetry {9 = 0) 
are just coordinate singularities, as is well-known, we can avoid them by chosing other 
coordinate systems. 

This metric has an Abelian two-parameter group of isometrics. We can choose two 
independent KVFs in the following way 

^ d ^ ( 2Mr\ , 2Marsm^9 , , , 

d , 2Marsin'^9 , sin^ 6* , 

^. = ^' = J. dt + ^^d^, (49) 

where = (r^ + a^Y — a^Asin^ 9. Now, we will study the Papapetrou field associated 
with the timelike KVF, ^t, 
_ 2M 



r - a cos 9) [dt - a sin 9d(fi A dr 



2Marsm(29) ^ / , . o , \ 

+ T^^de A [adt - (r^ + a^)difi] . (50) 

It is important to point out (see @]) that this Papapetrou field is proportional to the 
electromagnetic field of the Kerr-Newman metric (see also Q for a detailed study of 
this electromagnetic field). As we can see from (|48|) , the norm of is given by 

J, . (51) 

As is clear, it vanishes on the surfaces of infinite redshift (A = a? siv? 9), and in principle 
we cannot apply our procedure there. However, this situation is different from that 
of [Appendix A| , and later we will see that it is possible to use the general procedure of 
Sec. |. 

In the regions where the BL coordinates are valid, we find the following expressions 
for and a; 

^ = ^ |(r2 - a2 cos2 _ ^2 sin(20)rd6/} , 

2Masin6' r ^ , , n n o ,^-1 
uj = jcot 9dr + (r^ - cos^ 9)dej , 

and then, we can write the Ernst potential as follows 

2M, , . 2M 
= — (r + ta cos ( 



p2 r — ia cos 9 
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The eigenvalues of the Papapetrou field (|50| ) are given by 

2M S"^ 

a + i6 = S-, : rrrr = s— — , (52) 

^ (r-za COS 6')2 2M ^ ^ 

and therefore, it is regular everywhere excepting in the ring singularity = 0, where it 
diverges. On the other hand, it is remarkable that a + if3 is an analytic function of the 
Ernst potential. As we have pointed out before, this can only happens when dE = 
[see equation (0)], like in this example, where the 1-form E is given by 

E = s{dr + ia sin 9de) =^ dE = . 

Moreover, 

En = sdr, E'i,E^ = ^, Ej = sasinOde , E'}Ej^ = ^ , 

and hence, excepting in r = r± (A = 0) and in the axis of symmetry [9 = 0), the 
Papapetrou field corresponds to the generic case (Subsection p.l.l|) . In order to compute 
the principal directions of the Papapetrou field (|50|) in the generic case, we need to 
compute the 1-form D 

aAsin2 6' , ^ a^A sin^ 6'(A - sin^ 6*) 

p" 

Then, P and Q are 

P = (dt -asin^edcp) Q = —iadt - (r^ + a^)dcp) , 

p2 p2 



and finally, the principal null directions of the Papapetrou field are [see (|33|) ] 

J asin^M^ ^ A ^ asin^^A^ , ^, 

K = -dt + sdr H dcp , X = -dt — sdr H dcp , (53) 

p2 p2 pi pi 

which are just the multiple principal directions^]] of the Kerr space-time (see for 
instance |P), whose Petrov type is D. We can check, using Theorems 1 and 2, that 
K and L are indeed geodesic and shearfree. 

On the other hand, in [Appendix B| we have given the connection with the eigenray 
formalism of Perjes, showing that the existence of a geodesic (and shear-free) principal 
direction of the Papapetrou field is equivalent to the existence of a geodesic (and 
shear-free) eigenray. This result, together with the fact that the Kerr space-time is 
the only asymptotically-fiat stationary vacuum solution with geodesic and shear-free 
eigenrays [0, shows that the Kerr metric is the only one with such characteristics in 
which the alignment of the principal directions of the Papapetrou and gravitational 
fields takes place, which has been shown recently by Mars |jl^. In this paper, the 
characterization is given in the whole Kerr manifold, whereas the result of Perjes [^ 



was found only in the region where the KVF is timelike. However, taking into account 
that the eigenray formalism can be used also in the case of spacelike KVFs, Perjes' 
result could be extended to the whole manifold. 

% Note that the arbitrary sign s only serves to change from K to L and the converse. 
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Now, in order to show the capabihties of our formahsm, we will extend the analysis 
to the symmetry axis, the even horizons, and the surfaces of infinite redshift. On the 
axis of symmetry (where Ej = and Ef^Ej^^ ^ 0) we can find the principal directions 
by using expressions (|3^). However, the BL coordinates are not good on the axis. In 
order to solve this problem we can make the calculation in Kerr-Schild coordinates 
(see 0)- It is easy to see that the principal directions obtained can be expressed in 
terms of the coordinates {t,r}, they are given by 

K = -il- —— )dt + sdr, L = -[l- —— dt - sdr . 

The surfaces of infinite redshift (A = a? sin^ 6) constitute a controversial point in 
this example, because there the norm of vanishes [see (51)]. However, a + z/3 7^ 
(0) is different from zero on these surfaces, which tells us that we are not in the case 
of [Appendix A| . Then, we can follow the general procedure explained in Sec. 0. In this 
case, the Papapetrou field belongs to the generic case, and we can see that the principal 
directions are given by (|53|) . 

In the case of the event horizons, we are going to consider the maximal extension 
of the Kerr metric. As is well-known, it can be constructed by using advanced and 
retarded Kerr coordinates {V±, r, 6*, 99^} (—00 < Kt < 00, —00 <r<oo, 0<^<7r, 
< < 27r) 

dV± = dt± (r^ + a^) A^Mr , dif± = dip ± aA~Mr . 

Now, the line element is analytic (apart from regions I, II, and HI) also in r = r±. 
Therefore, we can already compute the principal directions on r±, where we have 
^R^Ra — 7^ 0- Taking into account the procedure explained in Sec. 0, 

one of the principal directions is given by En = dr, and the other one is obtained by 
contracting (^ with any timelike vector field. In our case, we find that on the event 
horizons ^ is given by 

:F = 2{dV± - a sin^ 6'dcp±) O dr , 

and therefore, the other principal direction is given by dV± — a sin^ 9d(p± (the vector 
field associated with this 1-form is d/dr, which is tangent to the curves V± = constant). 

To sum up, the maximal extension of the Kerr metric (see |TB[ for details), in the 
case < M^, is made up of an infinite chain of asymptotically-fiat type I regions, 
connected to type II regions (which contain trapped surfaces) and asymptotically-flat 
type HI regions (which contain the ring singularity) by means of the null hypersurfaces 
r = r±. We have just seen that we can compute the principal directions of the 
Papapetrou field (^) for the whole Kerr manifold, and in addition, that they coincide 
with the repeated principal directions of the Weyl tensor. For the sake of completeness, 
we want to remark that these results can be trivially extended to the cases = 
(where r+ = r_ and there is no region II) and > (where A > everywhere). 

We finish this example with a comment on the other KVF, (^9]). Both KVFs are 
privileged since is the only one (up to a multiplicative constant) which is timelike at 
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arbitrary large positive and negative values of r, and is the only one which vanishes 
on the axis of symmetry [9 = 0) and satisfies the regularity condition there (see j^). 
However, we can check that the principal directions of the Papapetrou field associated 
with do not coincide with the principal directions of the Weyl tensor. In this sense, 
our study privileges the KVF over the KVF ^<^. 

The next example will be an algebraically general vacuum space-time. One of the 
most simple families of this type is the following class of Kasner metrics (see also pf) 

ds^ = -dt^ + t^P^dx^ + t^P'^dy^ + fP'dz^ , 

where pi, p2 and are constants satisfying the following relationships 

Pi + P2 + P3 = 1 , pI+pI+pI = ^- 

This metric has an Abelian G3 group of isometrics acting on the spacelike hypersurfaces 
{t = constant}. We can choose one of these KVFs, for instance 

L = ^, ^cc = t^^'dx =^ F = 2pifP'dt A dx . 
ox 

Then, 

N = t^pi , ^ = -2pit^P'-^dt , w = , 

which means that the Ernst potential associated with $,x is simply S = —t^P'^ . Now, if 
we compute the eigenvalues of the Papapetrou field 

a = 2spitP^-^ , /3 = , 

we realize that it is a regular electromagnetic field, but it does not belong to the generic 
class, since 

E = -stP^dt Er = -stP'dt , Ei = 0. 

That is to say, we are in the particular case Ej = 0, and therefore the principal directions 
of the Papapetrou field are given by 



^ 9 d ^ d d 

The complex vector m which determines the orthogonal 2-planes to the principal 
directions, which is obtained from (P) or (|35|) , can be taken as follows 



4f 



9 . -,_„. 9 \ 



Now, we can check that the null principal directions (|54D of the Papapetrou field satisfy 
the condition of the Theorem 1, and therefore they are geodesic. In the general case, 
they do not coincide with the Weyl principal directions. Moreover, they satisfy the 
condition of the Theorem 2 only when P2 = Ps, in which case they are also shearfree 



and hence, via the Goldberg-Sachs theorem |16], the space-time is algebraically special 
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(specifically, Petrov type D), being ( [5^ ) the repeated null principal directions. Finally, 
using expressions (|i6|), we find that the expansion of k and £ is given by 

(the signs + and — are for k and £ respectively), and the rotation is zero since the 
rotation of $,x also vanishes [see (|15|)]. 

Another interesting example is a subfamily of the plane-fronted gravitational waves, 
pp waves (see for instance [§]). The line-element of these space-times is 

ds^ = -2dudv + 2dC,dC - 2Hdu^ , 

where 

that is, this is the subclass of the vacuum plane waves with the additional KVF 

L = -^, ^^=-dv- 2Hdu , N = -2H = -2 (/(C) + /(C)) • 

(These metrics have always d/dv as a KVF but it is a null KVF; see Appendix A). 
The Papapetrou field associated with is given by 

F = 2duA [H(dC + H^dC) , (55) 
Then, the 1-forms i/? and are given by 

ip = 2 [H(dC + H^dC) , u; = -2i [Hc^C - ^c^C) • (56) 

From these expressions we can write the Ernst potential associated with as follows 

S = 4/(0 . 
Furthermore, from (|56|) we deduce that 

Then, the Papapetrou field (|55| ) is singular {a + i(3 = 0), and it corresponds to the 
subcase (i) of the subsection Using (PI), the principal null direction is given by 

d 

k = IQHHcH?— , (57) 

which is the propagating direction of the plane wave, orthogonal to the wave fronts. 
Moreover, it is proportional to the null KVF d/dv and hence, taking into account the 
Theorem 3, k is geodesic and shear- free. From the Goldberg- Sachs theorem |]16|, we 



deduce that it is also a principal null direction of the space-time. In this case the metric 
is of the Petrov type N, with only one multiple principal direction, which is given by 
(^). Finally, the complex divergence p vanishes since d/dv is a constant vector field. 
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6. Remarks and conclusions 



In this paper we have developed a formalism for the study of the Papapetrou fields 
in vacuum space-times. In the first part of this paper, we have determined the null 
principal direction(s). In the case of regular Papapetrou fields, this is very useful in 
order to solve Maxwell's equations in curved space-times, since when we write them in a 
Newman-Penrose basis adapted to these principal directions the conditional differential 
system for them is second-order, contrary to what happens when we write them in an 
arbitrary Newman-Penrose basis, in which case it is a third-order differential system. 
In the second part of this work, we have studied the main differential properties of 
the null principal direction(s) of the Papapetrou field. Specifically, we have found the 
condition that a principal direction must satisfy in order to be geodesic. Moreover, for 
geodesic principal directions we have found explicit expressions for the optical scalars. 
This study allows us to study when a principal direction of the Papapetrou field is also a 
principal direction of an algebraically special vacuum space-time. Furthermore, taking 
into account the simplicity of the expressions for the differential properties, and the 
fact that they only depend on the derivatives of the Ernst potential and the principal 
direction, they provide a interesting way of introducing new Ansatze for the search of 
exact solutions. The examples studied here can serve as a guide. As we have seen, the 
principal directions of the Kerr space-time are aligned with those of Papapetrou field 
associated with the timelike KVF (see also |TB|, |T1|), which gives another characterization 
of this metric |T^, and the other two cases (Kasner and pp waves) are also examples of 
how the algebraic structure of the Papapetrou field can be adapted (for some KVFs) to 
the algebraic structure of the gravitational field. In this sense, the expressions for the 
expansion and rotation scalars of the principal directions are also useful. 

Finally, this study can be extended in several ways. One of them is to study 
the relationship between the principal direction(s) of the Papapetrou fields and the 
possible algebraic types of the space-time. Here, there are several cases depending on 
the multiplicity and degree of alignment of the principal direction(s) of the Papapetrou 
field On the other hand, taking into account that most of the properties studied 



can be expressed in terms of objects related directly to the Ernst potential, it seems 
reasonable to extend this study to other space-times in which the Ernst potential can be 
defined, for instance to Einstein-Maxwell space-times possessing a Killing vector field. 
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Appendix A. The case of a null Killing vector field 

In this Appendix we sum up briefly what happens when we consider a null KVF. In order 
to study this case, it is enough to remember that all the vacuum space-times admitting 
a null Killing vector fleld are known, they can be found in ^ (section 21.4; see also 
Ref. The Killing equations (p!oD tell us that the KVF is a geodesic, shear- free, and 



expansion-free null vector fleld, and the Einstein fleld equations imply the vanishing of 
the twist $, A d$, = 0. As it was shown by Dautcourt ||26|, these space-times can be 
divided into two different classes: 

-Class I is the family of the pp waves (see P]), in which the KVF is a constant null 
vector fleld {^a;b = 0), and therefore the associated Papapetrou fleld is identically zero. 

-Class II is determined by the following line element (see |^) 

ds"^ = -2xdu{dv + Mdu) + x~^/'^{dx^ + dy^) , 
where the function M = M{u,x,y) satisfles the following partial differential equation: 

(xM^),^ + xMyj, = 0. 
In this case, the null KVF is 

1=^, ^ = -xdu. (Al) 

Then, the associated Papapetrou fleld is simply F = du A dx. As is clear, it is singular 
and the corresponding principal direction is given by the null KVF (|A.1| ). Therefore, 
the principal null direction is geodesic and shear-free, and if we choose the KVF ( |A.1| ) 



to be the tangent vector fleld, the integral curves are affinely parametrized, being v an 
affine parameter. 

Appendix B. Connection with the eigenray formalism 



In reference Perjes developed a spinor calculus for stationary space-times from 



which, a triad formalism similar to the Newman-Penrose formalism was put forward. 
An applications of this technique is the search of exact solutions of Einstein's equations. 
To that end, Perjes introduced the concept of eigenray. Here, we show the connection 
between this concept and our development. 

Let us consider a non-null KVF ^ and a null vector fleld k normalized by k"'^a = 1- 
These two objects determine a unit vector fleld n by 

n=^\{k~N-'^), n'^n, = -sgn(A) , (B.l) 

where |x| and sgn(2;) denote the absolute value and the sign of x respectively. Now, we 



can show |]15[ that n is tangent to geodesies in the 3-space of the Killing orbits, which 
are called eigenrays, when it satisfles the following algebraic condition 

p'^Vfe + e\cn'u' = , Babe = mr'^'VdabcC'' , (B.2) 

where Pab = Qab + sgn(A)nanf, is the orthogonal projector to n. It can be also proven 
that condition ( [B.2|) is equivalent to say that k is geodesic (see [^), therefore n is 
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geodesic if and only if k is geodesic. Moreover, Perjes showed that a geodesic eigenray 
is also shear-free if and only if it is the null vector field associated through ( |B.1| ). 

Now, let us consider the Papapetrou field associated with ^. It is possible to show 
that the null vector field k in ( |B.1| ) is a principal direction of the Papapetrou field if 
and only if condition (|B.2[ ) holds. Therefore, this result shows that the existence of a 
geodesic eigenray is equivalent to the existence of a geodesic principal direction of the 
Papapetrou field (and the same for the shear- free case). 
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